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An avalanche photodiode (APD) array for ground-based optical communications 
receivers is investigated for the reception of optical signals through the turbulent 
atmosphere. Kolmogorov phase screen simulations are used to generate realistic 
spatial distributions of the received optical field. It is shown that use of an APD 
array for pulse-position modulation detection can improve performance by up to 
4 dB over single APD detection in the presence of turbulence , but that photon- 
counting detector arrays yield even greater gains. 


I. Introduction 

Ground-based reception of optical signals from space suffers from degradation of the optical phase 
front caused by atmospheric turbulence, leading to a reduction in the effective diameter of the receiving 
telescope and to random fluctuations of the point-spread function in the focal plane. A proportional 
increase in the receiver’s field of view in order to collect all of the signal also causes a corresponding increase 
in the amount of interfering background radiation, resulting in degraded communications performance. 
These problems may be mitigated through the use of an optical detector array assembly in the focal 
plane that can adaptively select areas of higher signal density while ignoring areas predominated by 
background noise. In [1], the performance of a detector array composed of photon-counting detectors 
was evaluated and found to yield up to a 5 dB improvement over a conventional single-detector photon- 
counting receiver. However, the current baseline receiver design for deep-space optical communication 
utilizes readily available avalanche photodiode detectors (APDs) rather than photon-counting detector 
arrays, which are still in the development stage at the wavelengths of interest for optical communications. 
In this article, we extend some of the results obtained in [1] to the APD array case. 


II. APD Output Model 

In [1], the derivations of signal models for the single photon counter and the photon-counting array were 
given in considerable detail in order to obtain the optimally weighted photon-counting array performance. 
It was then shown that a suboptimal array (referred to as the adaptive synthesized detector or 0-1 
subarray) consisting of the optimal number of unweighted array elements (i.e., weights of zero or one) 
yielded performance very near to that of the optimally weighted array. In this article, we present results 
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for the 0-1 APD subarray. In our discussion of APD output statistic modeling, we follow the formulation 
found in [2]. 

The average number of photons absorbed by an APD illuminated with total optical intensity A (t) in 
T s seconds can be expressed as 


n = [ A (t)dt (1) 

hv J o 

where h is Planck’s constant, v is the optical frequency, and r] is the detector’s quantum efficiency, 
defined as the ratio of absorbed to incident photons. The actual number of photons absorbed, n, is a 
Poisson-distributed random variable with probability 


P(n\n) = ^ e- (2) 
n\ 

The probability p(m\n) that an APD generates m output electrons given n mean absorbed photons can 
be derived from the Mclntyre-Conradi distribution [3], but may be approximated by the continuous 
Webb density function [4] or by a Gaussian density [2,5]. Added to the random number of APD output 
electrons is an independent Gaussian thermal noise charge from the follow-on electronics [2]. For M-ary 
pulse-position modulation (PPM) with slot duration T s , the total charge is integrated over each slot time, 
T s , resulting in a vector of M independent observables for each received PPM word. It was shown in [6] 
that, given these observables, the maximum-likelihood detector structure consists of choosing the PPM 
symbol corresponding to the slot with the largest accumulated charge value. 

In this article, we use the Gaussian approximation for p(m\n) in order to simplify calculation of PPM 
symbol-error probabilities. This approximation is adequate for high background conditions, as shown in 
[5]. The slot statistic consisting of the sum of APD output electrons and amplifier thermal noise is thus 
also Gaussian, leading to M-ary PPM symbol-error probability given by 
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where a 2 ) is the Gaussian density function with mean fi and variance cr 2 , and 4>(x) is the 

Gaussian distribution function. The background and signal means and variances are given by pb = 
qGn b + / S T S , n s = qG(n b + n 3 ) + / S T S , = [2 q 2 G 2 Fn b + 2 qI s T s + (4 k,TT 8 /R l )\ BT s , and cr 2 = 
\2q 2 G 2 F{fib + n s ) + 2 qI s T s + (4 kTT s /R)\ BT s [2]. Here, F = kG + (2 — 1/G)(1 — fc), G is the aver- 
age APD gain, k is the APD ionization ratio, q is the electron charge, k is Boltzmann’s constant, T is 
the equivalent noise temperature, B is the single-sided noise bandwidth, R is the load resistance, and I s 
is the APD surface-leakage current. The mean numbers of absorbed background and signal photons are 
given by fib and n s , respectively. The noise bandwidth, H, is set to that of an ideal T s second integrator, 
i.e., B = 1/2 T s . The APD bulk dark current is modeled as part of the background radiation, i.e., it is 
incorporated into the value of 


III. Adaptive Synthesized Array 

If we consider a rectangular array of detectors consisting of K x L detector elements, the optical signal 
intensity incident upon each array element may be denoted by A where 1 < i < K and 1 < j < L. 
The average number of photons incident upon the ijth detector is then fiij = (rj/hu) f^ s A ij(t)dt. Using 
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the APD modeling described above and assuming that each array element observes the sum of a signal 
field plus multimode Gaussian noise field with an average noise count per mode much less than one, the 
outputs of the APD array may be modeled as conditionally independent Gaussian processes, conditioned 
on the average signal intensity over each detector element [7,8]. Given the APD array-element observables, 
the optimum maximum-likelihood detection scheme may be derived, consisting of a weighted sum of array 
outputs [1]. We consider a simpler real-time suboptimum detector whose array weights are either zero or 
one, i.e. , the array-element outputs are either included or excluded at any given time in making a PPM 
symbol decision. We therefore list the detector elements in decreasing order of average signal intensity, , 
and compute the probability of error for the first detector element’s signal intensity plus the background 
incident upon that element, then form the sum of signal energies from the first two detector elements 
(plus background for two detector elements), and so on, until the minimum error probability is reached. 
Each set of detectors may be effectively considered to be a single detector, so that no weighting is applied 
to account for variations in the signal distribution over the detector elements included in that set. The 
set of detector elements that achieves the minimum probability of error is the best “synthesized single 
detector” matched to the signal-intensity distribution. As in [1], a sample signal-intensity distribution may 
be generated using Kolmogorov phase-screen algorithms [9]. Note also that the process of optimizing the 
subarray requires calculation of the PPM symbol-error probability for each number of detectors, as well 
as knowledge of signal and background energy levels. A practical real-time implementation of this array 
detection scheme will involve additional parameter-estimation algorithms that have yet to be developed. 

In evaluating the performance of an APD array, we must use an appropriate model for the thermal-noise 
contribution due to follow-on electronics in order to arrive at appropriate values for the load resistance, 
i?, constrained by maintaining sufficient bandwidth to reproduce the signal pulse. We make the following 
assumption: the detector elements are followed by individual transimpedance amplifiers, each of which has 
a feedback resistance whose value is determined by R = 7 x 10 12 T S ohms, as per the FOCAS optical link 
budget program [10]. The outputs of these amplifiers may then be sampled and added together digitally 
to obtain the desired signal, as illustrated in Fig. 1. With fine enough quantization, the sampling and 
summing operations may be considered to be essentially noiseless for this application. Note that the 
total variance due to thermal noise at the output is the sum of the individual thermal-noise variance 
contributions from the load (i.e., feedback) resistors following each of the detectors. 

Assuming the above model, the charge accumulated during each time slot for the ij th detector element 
is Gaussian distributed with mean fiij = qGfiij + I S T S and variance of- = q 2 GFfiij + qI s T s + 2kTT s / R j, 
where is given by n s ^j + fi b for a signal slot and n b for a noise slot. Since the outputs from the detec- 
tor elements are independent Gaussian random variables, the composite signal from adding N detector 
outputs belonging to subset S has mean 
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where S denotes the set of detector elements selected for best performance, and where we have as- 
sumed that each transimpedance amplifier has the same feedback resistance value of R = 218, 750 ohms 
as determined from FOCAS for the required signal bandwidth of 32 MHz. Therefore, as the number of 
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Fig. 1. Load resistance models for an APD array (where 
resistors are denoted by "R," amplifiers by "A," and ana- 
log-to-digital converters by "A/D"). 


detectors, TV, in the subarray increases, the effective total load resistance used for computing the variance 
of the combined output varies as R/N. The other parameters used in the calculations were G = 100, 
k = 0.07, T = 300 K, T s = 31.25 ns, r] = 0.4, and I s = 3.9 x 10" 10 nA. 

IV. Numerical Results 

In order to generate a sample function of the spatial distribution of the signal incident upon the detector 
plane, sample fields were generated using a Kolmogorov phase-screen program [9] with an atmospheric 
correlation length of ro = 4 cm, representing moderate turbulence, resulting in a matrix of complex signal 
amplitudes. The field intensity generated in the detector plane by the simulation then was integrated over 
the elements of a 16 x 16 = 256 detector array that was dimensioned to encompass the significant extent 
of the signal distribution in the detector plane, corresponding roughly to an optimized single-detector 
design. A constant average background photon energy of n 5 is assumed over each detector element. For 
a given sample function of the intensity distribution, the 256 detector elements were sorted in decreasing 
order of average signal energy, and M-ary PPM symbol-error probabilities were calculated for increasing 
numbers of detectors, starting with the detector with the highest incident signal energy [1], in order to 
specify the adaptive synthesized detector. Note that the results presented here are for a single realization 
of the turbulent signal distribution. However, in [1] it was observed that different turbulence realizations 
(with the same atmospheric correlation-length parameter value) did not result in significantly different 
array detection results for the photon-counting channel. Because the photon arrival process is the same 
for both models, it is reasonable to assume that this will hold for the APD detector array as well. 

Performance results were obtained for two types of receivers: (1) the adaptive synthesized detector 
(the optimum number of unweighted detector elements are used) and (2) the single large detector (all 
256 detector elements are given unity weight, effectively synthesizing a nonadaptive single large detector). 
The results of the APD array simulations are shown in Figs. 2 through 9 and are compared with analogous 



Fig. 2. PPM symbol-error probabilities for 2-PPM, n b = 0.1 absorbed 
background photon per detector element per slot. 



Fig. 3. PPM symbol-error probabilities for 2-PPM, n b = 1.0 absorbed 
background photon per detector element per slot. 


photon-counting detector results. The APD results were obtained using a Gaussian approximation for 
the distribution of the detector output, composed of Webb-distributed secondary electrons plus additive 
Gaussian noise. The photon-counting results assumed Poisson-distributed observables, which accurately 
represent the statistics of multimode signal plus background radiation under typical operating conditions 
[8]. PPM modulation with M — 2, 4, 16, and 256 was investigated, operating in the presence of low-to- 
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Fig. 4. PPM symbol-error probabilities for 4-PPM, n b = 0.1 absorbed 
background photon per detector element per slot. 



Fig. 5. PPM symbol-error probabilities for 4-PPM, n b = 1.0 absorbed 
background photon per detector element per slot. 


moderate background levels that generate an average of 0.1 absorbed photon per detector element over 
each slot time (or an average of 25.6 absorbed photons per slot over the entire array, which can also be 
viewed as a single large detector), as well as high background levels of 1 absorbed photon per detector 
element per slot (or 256 absorbed photons per slot over the entire array). For the 256-PPM photon- 
counting results shown in Figs. 8 and 9, accurate computation of the exact symbol probability was 
prohibitively complicated, so a union bound (as discussed in [1]) was used. 
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Fig. 6. PPM symbol-error probabilities for 16-PPM, n b = 0.1 absorbed 
background photon per detector element per slot. 



Fig. 7. PPM symbol-error probabilities for 16-PPM, n b = 1.0 absorbed 
background photon per detector element per slot. 


The results for low-to- moderate background levels (0.1 absorbed photon per detector element per slot) 
are shown in Figs. 2, 4, 6, and 8 for each value of M. The gain from APD array processing over a single 
APD is about 3.6 dB at a symbol-error probability of 0.001 for all values of M. It is apparent that for the 
case of low-to- mo derate background levels further improvements are possible with the use of an adaptive 
APD array to reduce the size of the subarray, thus limiting the background contribution to the total noise 
process. 
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Fig. 8. PPM symbol-error probabilities for 256-PPM, n b = 0.1 
absorbed background photon per detector element per slot. 



Fig. 9. PPM symbol-error probabilities for 256-PPM, =1.0 

absorbed background photon per detector element per slot. 


When operating in high background environments characterized by an average of 1 absorbed photon 
per detector element per slot (or 256 absorbed photons for the entire array, or for a single large detector), 
the improvements due to APD array processing tend to be greater. At a symbol-error probability of 
0.001, adaptive APD array processing now yields approximately 4 dB over the single large APD for all 
cases considered. The reason is that the background contribution to the total noise variance in Eq. (5) 
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dominates; reducing the sum of the elemental background contributions by restricting the observation 
region to a much smaller subarray therefore results in a proportionally smaller noise variance. 

While significant, the improvements demonstrated by the use of adaptive APD arrays over a single 
large APD detector are nevertheless smaller than gains achievable through the use of adaptive photon- 
counting detector arrays [1], as shown in Figs. 2 through 9. The reason is that photon-counting detectors 
are not hampered by excess noise due to the avalanche multiplication process or by thermally generated 
circuit noise, which tend to limit the performance of conventional APD detectors. For the case of low-to- 
moderate background levels, the use of adaptive photon-counting arrays results in gains of 7 dB to 7.5 dB 
over the APD array, or 10.6 dB to 11.1 dB over the single large APD, at symbol-error probabilities of 
0.001. For the high background case, the corresponding improvements are 5 dB to 5.5 dB over the APD 
array and 9 dB to 9.5 dB over the single large APD. Adaptive photon-counting arrays therefore provide 
greatly enhanced capabilities for the reception of PPM-modulated optical signals through the terrestrial 
atmosphere. 


V. Summary and Conclusions 

The results of earlier studies concerning the use of optical focal-plane photon-counting detector arrays 
together with adaptive signal processing to combat the effects of atmospheric turbulence [1] have been 
extended to include APD arrays used in a similar configuration. The total detector output process, 
composed of APD secondary electrons due to absorbed photons together with additive circuit noise of 
thermal origin, was modeled as a Gaussian random process with mean and variance consistent with the 
statistics of the component processes. Adaptive array performance was evaluated for operation in both 
low and high background environments for the case of PPM-modulated received fields, and compared 
with the performance of the photon-counting adaptive array previously documented in [1]. It was found 
that while adaptive APD detector arrays provided improvements of 3.6 dB to 4 dB over the performance 
of single large APD detectors designed to collect most of the received signal, the use of photon-counting 
arrays instead of single large APD detectors provided much greater additional improvements, amounting 
to as much as 11 dB in low-to-moderate background environments and up to 9.5 dB in high background 
environments. Improvements of this magnitude warrant serious consideration of the adaptive photon- 
counting array technology for future DSN optical communications applications. 
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